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Nonlinear light propagation in a single-mode micron-size waveguide made of semiconducting exci- 
tonic material has been theoretically studied in terms of exciton-polaritons by using an analysis based 
on macroscopic fields. When a light pulse is spectrally centered in the vicinity of the ground-state 
Wannier exciton resonance, it interacts with the medium nonlinearly. This optical cubic nonlinear- 
ity is caused by the repulsive exciton-exciton interactions in the semiconductor, and at resonance 
it is orders of magnitude larger than the Kerr nonlinearity (e.g., in silica). We demonstrate that a 
very strong and unconventional modulational instability takes place, which has not been previously 
reported. After reducing the problem to a single nonlinear Schrodinger-like equation, we also ex- 
plore the formation of solitary waves both inside and outside the polaritonic gap and find evidence 
of spectral broadening. A realistic physical model of the excitonic waveguide structure is proposed. 
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I. INTRODUCTION 

The nonlinear phenomena of modulational instability 
(MI) and solitons formation make the basis of modern 
nonlinear fiber optics, where the Kerr (cubic) nonlinear- 
ities are mainly used.'^' MI consists in the exponential 
("parametric") growth of equally spaced spectral side- 
bands from a continuous wave pump, and is equivalent 
to the degenerate four- wave mixing.^ It is a particularly 
useful phenomenon in microstructured fibers, where it 
allows the generation of coherent light at frequencies 
where lasers are not easily available.El The MI instabil- 
ity is strictly connected to the existence of solitons, or 
solitary waves, i.e. localized short pulses that are invari- 
ant in the propagation, and the shape of which is de- 
termined by the spe cific nonlinearity and the dispersive 
features of fiber.lEI Plane waves tend to become unsta- 
ble due to the fact that the actual "normal modes" of 
the fiber are nonlinear optical solitons.^i^ A particularly 
useful phenomenon connected to solitons is the so-called 
supercontinuum generation, which is the sudden and ex- 
plosive spectral broadening of an input energetic short 
pulse in the fiber.l^ 

However, optical fibers are not the only media where 
MI and solitons can be observed and used. The gen- 
eral investigations of optical nonlinear effects in di- 
electrics and semiconductors also have history of sev- 
eral decades,!^^ and recently many applied researches 
have been started to investigate the potential of the soli- 
tary wave sector of semiconductor materials near exci- 
tonic resonances, for example, in semiconductor micro- 
cavities or semiconductor gratings. The dominant 
optical cubic nonlinearity for the polarization field in ex- 
citonic materials, which we investigate here, is caused by 



the exciton-exciton Coulomb repulsive interactions and 
essentially differs from the cubic Kerr nonlinearity for 
the electric field (e.g., in silica). The excitonic nonlin- 
earity is orders of magnitude larger than the Kerr one, 
but it also possesses a resonant nature - as we will see in 
the following, it substantially weakens when the detun- 
ing of the central frequency of the pulse from the exciton 
resonance increases. The first difference allows one to ex- 
pect the observation of MI and solitons formation not on 
the kilometer (e.g., like in silica fibers) but on microm- 
eter scale and for lower pulse intensities. This opens up 
the possibility of on-chip integration of optical devices 
for nonlinear frequency conversion based on waveguides 
made of semiconducting excitonic materials. 

In this work we theoretically study all the mentioned 
optical nonlinear properties of such excitonic waveg- 
uides. In Sec. In] on the basis of known nonlinear dif- 
ferential equation for the complex envelope of macro- 
scopic medium polarization fielcP and Maxwell's equa- 
tions for the light pulse we derived a general set of equa- 
tions, which govern the light propagation in the excitonic 
waveguide. Further, we applied this set of equations 
to investigate the instability of long light pulses prop- 
agating in the excitonic waveguide with respect to small 
perturbations and without any approximations discov- 
ered a very strong modulational instability of the pulse 
spectrum with enormous gain (Sec. III). This is com- 



pletely confirmed by the direct numerical simulation of 
the pulse spectrum evolution using the general set of 
equations. Next, we reduced the general set of equa- 
tions to a single nonlinear differential equation for the 
polarization field envelope, which in its turn under the 
slowly varying envelope approximation (SVEA) was re- 
duced to the nonlinear Schrodinger equation with an ex- 
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act and complete set of coefBcients. Within the so-called 
"tail analysis" formalisnlSi all the parameters of corre- 
sponding solitonic solution of this equation have been 
obtained for light pulses spectrally centered both outside 
and even inside the polaritonic gap, where the propaga- 
tion of linear plane waves is forbidden, while solitons can 
nevertheless form and propagate if the incident pulse in- 
tensity is large enough (Sec. IV). In this way we have 
substantially generalized previous theories in the exist- 



ing literatures^ without any additional approximations. 
These analytical derivations are also supported by direct 
numerical simulations using the general set of equations 
under SVEA. Finally, the obtained results allowed us to 
propose a realistic physical model of the excitonic waveg- 
uide and discuss the experimental conditions for the ob- 
servation of above-mentioned effects using GaAs as the 
representative excitonic material (Sec. \V\. 



II. NONLINEAR EQUATION FOR 
POLARIZATION 



In Ref. [5] it was shown that the set of semicon- 
ductor Bloch equations for the excitonic medium in 
the low excitation regime can be reduced to a sin- 
gle nonlinear differential equation for a macroscopic 
complex envelope P{z,t) of the medium polariza- 
tion 7'(r,i) = {P{z, t)M{r j_)exp[ikz - iut] + c.c.} /2 
coupled to the external electric field f(r,<) = 
{E{z, t)F{r±) exp[ikz — iut] + c.c.} /2, where z is the 
longitudinal coordinate and rj^ are the transverse coor- 
dinates. Here we assume the existence of a fundamental 
guided mode, thus allowing a factorization of the dimen- 
sionless transverse modal distributions M and F from 
their respective envelopes. Maintaining the most essen- 
tial nonlinear terms, after the separation of variables this 
equation can be written in the following form; 



idtP + {Alo^ + i-f^)P ~ aTi\P\^P 

+ aeB^(rMF-/3r2|Pn =0, 



(1) 



where = uj — uj^ is the detuning of the pulse 

spectrum central frequency lu from the Is-exciton 
resonance Wx, 7x is the exciton damping parameter 
introduced here phenomenologically, a and /3 - non- 
linear coefficients, which are completely defined by 
the microscopic (quantum) properties of the exci- 
tonic medium (within the jellium two-band exciton 
model one has, as in Ref. HI a = 26u;b/3 and /3 = 7, 
where Wb is the exciton binding frequency), £b is 
the bulk background dielectric constant and d is the 
photon-exciton coupling parameter. It is well-known 
that the photon-exciton coupled state, the so-called 
polariton, shows a region of forbidden frequencies in its 
dispersion relation - the polaritonic gap, which width 
is given by d. Parameters F depend on the relation 



between F and M:_ Fi = J \M\^dr±/ J \M\^drj_, 
F2 = /|M|2MFdri//|M|2dri, Fmf = 
J MFdrj_/ J \M\'^drj_. The modal distributions F 
and M are localized functions,'^ therefore Fi < 1, 
F2 ^ 1, and Fmf ^ 1, while Ffm ^ 1 as the polarization 
field modal distribution M is completely determined by 
that of the electric field F. The cubic term in Eq. ([T]) 
originates from the repulsive Coulomb exciton-exciton 
interactions in the semiconductor, while the square term 
is responsible for the phase-space filling effect (also 
called the Pauli blocking) when the density of excitons is 
comparable to the Mott density (for which ~ 1), 

what leads to a decoupling of the external electric field 
from the polarization field. 

In the low excitation regime considered in Ref. |8j for 
exciton densities far below the Mott density, one can ap- 
proximate Eq. ([1]) as: 

idtP + [Auj^ + i"f^)P ~ aFilPpP + dVuFSBE = 0. (2) 

We perform the following analytical derivations and 
numerical simulations on the basis of Eq. ([2) cou- 
pled to Maxwell's equations for the electric and mag- 
netic field envelopes of the incident pulse, 'H{r,t) = 
{H{z, t)F{r±) exp[zfc2; — iujt] + c.c.} /2. In the most gen- 
eral dimensionless form this set of equations is given by: 



8x7] = {iuj' - dT)i^J + FfmA((9) 
dxijj = {ii-^' — dT)Ti — ik'-tp, 



ik'rj, 

Dt^P = i ([i7x + - LP + FmfVj) 



(3) 



where the following redefinitions and scalings have been 
applied: ip — E/Eq, ip = P/Po, rj = H/Hq, x ~ z/zq, 
T = t/tQ, uj' = ujtQ, 7^ = 7x^0 , Awx = Awxto, A = ato, 
k' = fczo, and Pq = l/VlYato, = Po/{eBdto), Hq = 
nEo, zq — cto/n, c is the velocity of light in vacuum, n — 
y^EB is the non-resonant background refractive index, and 
to is an arbitrary time scaling parameter, which will be 
properly chosen in the following. 



III. MODULATIONAL INSTABILITY 
ANALYSIS 

We now use system ([s]) to analyze the linear stabil- 
ity of a long pulse (ideally a continuous wave) propa- 
gating in the excitonic optical waveguide with respect 
to small perturbations. This is the MI analysis, see 
also Ref. [H By substituting the perturbed field and 
polarization envelopes {il^;r];(p}{x,T) = ({V'o! ^70; ¥^0} + 
{a; g;p}{x, T)) exp[iqx] into system ([S]), setting the small 
perturbations as {a; g;p}{x,T) — {ai; gi;pi} exp[iKx — 
iST] + {02; 32; P2} exp[i(5r — inx] and imposing the solv- 
ability condition of system ([s]) one can obtain within the 
first order perturbation theory- the dispersion relation 
k((5) for perturbations: 
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FIG. 1. (a) Numerically obtained spectral evolution of a 
long pulse propagating along the a;-axis (according to (l3|, 
the initial pulse is 'i){x = 0,T) = Vo exp [-(2r/Tw)^"] , 
Tw = 50, m = 20, ipo = 10.2, and the other parameters 
are Au'^ = —1.1, 7^ — 0.11). (b) Corresponding analytical 
prediction on the MI gain G{5) peaks, given by Q. ui' is 
positioned at the origin. 



where k is the perturbation wavenumber, S is the detun- 
ing from the central pump frequency, and q = uj'[1 + 
ArMFrFM/(</5o ~ ^^DV^"^ ~ ^' ■ Those frequency regions 
where Im k((5) < correspond to the exponential growth 
of perturbations and this defines the MI gain spectrum 
G{5) — 2max |Im k(i5)|. Among all the solutions of 
Eq. Q, which is of the fourth order in k, we select, for 
a given 5, the one with the maximum absolute value of 
the gain. 

To confirm the above analytical results we performed 
direct numerical simulations of the propagation of long 
pulses in the excitonic waveguide by using system ([3| 
without any approximations and obtained the spectrum 
evolution (see Fig.[lja)) in complete accordance with the 
analytical predictions on the MI gain maxima positions 
(see Fig. [lib)). 

A fourth order Runge-Kutta algorithm for the mag- 
netic, electric and polarization fields has been employed 
with the following initial conditions: ^{x = 0,T) is a 
wide super-Gaussian pulse, 'q{x — 0,T) — ±^p{x — 0,T) 
(for respectively forward and backward propagation) and 
ip{x = 0, T = 0) = 0, which means that the medium is 
not polarized before the field arrival. We choose to = 1/a 
(for GaAs this is ~ 8 ps) as natural temporal scale of the 
system. Here we set k' = uj' and for the simulations 
removed the large number uj' ^ Lo^/a from ([3]) by mul- 
tiplication of the two first equations by A/w cind space 
rescaling zo zqX/lo' ^ c/{nuj) ^ 0.04 /^m, which is 
the general physical space scaling of the whole developed 
model (from now on GaAs is selected as a representa- 
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FIG. 2. Dependence of the MI gain G{5) on the frequency 
detuning Aio'-^ according to Q for the fixed value of external 
electric field amplitude tpo = ((^o/rMF)(v5o ~ ^'^x)- The used 
parameters are i/io = 5, 7^ ~ 0.019. 



five material, for numbers see Sec. |V). With this scal- 
ing it can be seen from Fig. [T] that the MI gain at its 
maxima reaches values ~ 100 /im~^, which is enormous 
in comparison with that for silica optical fibers possess- 
ing nonlinear Kerr effect with the MI gain ~ 10 km~^. 
This example clearly shows that the nonlinear phenom- 
ena, which manifest themselves, e.g., in optical fibers on 
the km-distances, can be observed in the excitonic waveg- 
uides when the incident pulse has propagated for just a 
few microns. 

It should be noted that in our model the excitonic 
damping 7^ has been naturally taken into account not 
only in the simulations but also in the analytical results 
by the redefinition Aoj^ — Aw^ + «7x- The increase of 7^ 
makes the MI peaks smaller in gain and spectrally wider. 
The dependence of MI gain on the detuning Aw^ with 
small fixed nonvanishing 7^ is very unconventional and is 
shown at Fig. [2j Peaks with unusually large gain appear 
in the spectrum because of the resonant nature of the ex- 
citonic nonlinearity, which, of course, affects the disper- 
sion relation (|4| - complex poles appear in its roots. The 
positions of these resonant gain peaks shown at Fig. [2] 
are completely determined by such poles, and the mag- 
nitude of the gain is consequently much larger than in 
non- resonant instabilities. In Fig. [2] also much lower cen- 
tral peaks can be seen, which are almost flat in the scale 
of Fig. [l]^b). These peaks shapes are analogous to those 
observed in the MI gain spectra of optical fibers in the 
vicinity of the pump frequency, where only non-resonant 
Kerr nonlinearity is present .'i' 
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IV. FORMATION OF SOLITARY WAVES 

We have also studied the sohtary wave formation under 
intense and short pulse excitation of the excitonic waveg- 
uide in the coherent regime (7x = 0). After the change 
of variables = x — VT (where V is the dimensionless 
soliton group velocity) system ^ in analogy with Ref.lTCl 
can be reduced to a single third order nonlinear differen- 
tial equation for the polarization envelope ip: 



(5) 



where: 



pi EE F(V^2 _ 

p-2, = ArV^ - (V^ - l)Aa;^ - 2V{VJ - fc'), 
P3 = 2Vuj'\v - 2I^uj'^{VJ - fc') + V^(fc'2 - 
P4 = (A;'2-c^'2)At^; + Arc^'2, 
CTi = _ 1^ = 2(ya;' - fc'), cts = fc'^ - cj'^, 

and Ar — AFMprFM- Under the slowly varying envelope 
approximatiorfS (SVEA), because p^+i ^ Pi, i — 1, 2, 3, 
and (Ti+i ^ (Ti, « = 1,2, Eq. ([5| can be reduced to the 



nonlinear Schrodinger equation: 

Pid\H^ - PW + <^iW?^ = 0, 



(6) 



where we also set p3 = as the definition of group ve- 
locity. As is well known,*^ Eq. (|6| has the solitonic 
solution: 



where d = ^Jlp/^jaj,, and h 
electric field takes the form: 



d sech , 

\/P2lP\ 



(7) 

Therefore, the 



= ^/Emf {-iVlh tanh - Ac^^, + . (8) 

It is important to note, that here SVEA has been used 
only after the reduction of the full system ([s]) to the sin- 
gle third order differential equation, that, unlike Ref. fTTl 
makes the coefficients in ([6]) exact and keeps the most 
essential cubic term at CT3. 

In order to explicitly obtain the dispersion relation and 
other parameters of the propagating pulse we apply the 
so-called "tail analysis" described, for example, in Ref. [51 
Starting from the full system ([s]) and looking for solutions 
in the form {V"; '7}(0 = {"00; 0o; exp[-^/l^], which 
reproduces the exponentially localized shape of a soliton 
for ^ — > 00 (at its "tail"), one can linearize the system ([3| 
and finally obtain: 



fc' = 



V = 



Cn 

cr 



^2_Yl_ (1 - Vl^')Aw; - 2{uj' - k'/V) 



_ i)A< + Ar 

Ar (cj'2 - 2w'Aa;; - l) 
^ 2w'(l + Aa;;2) 



d^ 

C2 



= 2 

_ 1 
~ 2 



AujL 



Ar 



e/eB - 1, 

4^-2 \ 1/2- 
C2+(Cl + ^ 



1 + Awi 







[0- 







(1 + A< _ Ar Ac^;;) 



2Ar 



(9) 



The set of Eqs. ^ completely defines all soliton parame- 
ters (the dispersion law k'{uj'), the group velocity V, the 
amplitude d, the spatial width parameter b) and contains 
just a single free parameter to, which in this Section is 
the soliton temporal width. All the other parameters are 
fully determined by the excitonic medium properties (cJx, 
a, Eb, etc.). 

Although the developed model is valid for relatively 
wide solitons {to ^ 1/a) , it g eneralizes previous theories 
present in the literature,'^'i^ and allows one to calculate 
all the soliton parameters for the incident pulses spec- 
trally centered both outside and even inside the polari- 
tonic gap, where solitons could also form if the incident 
pulse intensity is large enough (see Fig. [s]). Previously 
(e.g., in Refs.[6land[T3|). in- and out-of-gap solitons were 
treated separately with some additional approximations. 
It is clear from Fig. [3] that to form in-gap solitary waves 



one should input an intensity, which is orders of magni- 
tude higher than in the case of the out-of-gap solitons. 
The amplitude of the latter grows with the detuning be- 
cause far from the exciton resonance the nonlinearity is 
weaker (see Sec. |V]) and for the observation of nonlinear 
effects the input intensity should be increased. 

By using system ([s]) under the SVEA (this is analogous 
to what is done in Ref. [11] and valid only out of the gap, 
where k' — )■ ui' and (T2 tends to become the dominant 
nonlinear coefficient in Eq. ([s])): 

J dxip = jFFMA(y5/2, , , 

\ drP = t{H + ^< - 1^1'] ^ + TmF^) , ^ 

where t = T — x, we simulated the propagation dynam- 
ics of a short intense sech-shaped pulse in the excitonic 
medium and demonstrated solitary waves formation for 
the out-of-gap regime (see Fig. |4|^a)) and an essential 
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FIG. 3. (a) The amplitude of pulse soliton as a function 
of the detuning Aoj^ (in the units of dimensionless gap Ar) 
for the different soliton temporal widths to — s/a according 
to ([7]l-([9]|. The exciton resonance is positioned at the origin, 
the gap is indicated by the grey lines, (b) The correspond- 
ing temporal shapes of solitons both inside and outside the 
polaritonic gap in the case of s = 2. 
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FIG. 4. (a) The formation of well separated out-of-gap soli- 
tary waves is numerically demonstrated for the short intense 
sech-shaped pulse propagation in the excitonic medium, 
(b) The evolution of corresponding spectrum with the essen- 
tial broadening. Calculations are performed by using sys- 
tem ( |10[ ) and initial pulse — 0,t) = ipo sech(r), ipo = 20, 
while other parameters are = 1-2, 7^ — 0.19. 



broadening of the spectrum (see Fig. Qb)). In the case 
of w > Wx this spectral broadening, as well as the MI 
sidebands detuning from uj, should be limited by the ex- 
citon binding frequency value to avoid the excitation of 
band electrons. Here we employed a numerical method 
similar to that used in Sec. IIIIl 



V. REALISTIC EXCITONIC WAVEGUIDE 
PARAMETERS 

On the basis of above results we can discuss now the 
structure and parameters of a realistic excitonic waveg- 
uide and the physical conditions (and constraints) for the 
observation of considered effects by the example of some 
specific excitonic material. 

As it was mentioned previously, GaAs was selected 
as the representative semiconducting excitonic material, 
for which at low temperatures and for high-quality sam- 
pjj,g,i8 20 ^ i^i^ meV, /iWb = 4 meV, ha = 

0.08 meV, h-y^ ~ 0.03 meV, and £b = 12.56. The low 
temperature condition is menti oned here because, as it 
was shown expcrimcntallyf^^^ for GaAs the polaritonic 
effects take place only below ^ 20 K. However, for those 
samples and temperatures, for which ha < hjy^ (for ex- 



ample, in presence of inhomogeneous broadening), the 
polaritonic gap will not be visible and only out-of-gap 
effects can be observed. 

According to Refs. E] and [21] one can properly re- 
turn dimensionality to the macroscopic polarization 
field in the following way: P — >■ uP, where u = 
[2(icv/(7''aneo)][eoc«/2]^/^ (P is measured in the units of 
intensity^, dcv is the intcrband dipole matrix element, 
Oq is the Is-exciton Bohr radius, Eq is the vacuum permit- 
tivity, and the spin summation is also taken into account. 
Therefore, the nonlinear coefficients in Eq. ([T]) take form 
{a;/?} — {q;;/3}/u^, while the coupling parameter be- 
comes equal^ a = 2(i^y/(7raQfeoeB) in complete accor- 
dance with the polaritonic gap definition.^-' For GaAs, 
flo ~ 12 nm and dcv/e — 0.4 nm. These parameters 
allow one to calculate the corresponding nonlinear coeffi- 
cients: a ~ 2040 mV(W s) and ^ ~ 2.33 • 10"!" mVW. 
Now we can estimate the limiting field intensities from 
the condition \Pf - 1//3: P^^ ~ 0.43 MW/cm^ and 
El^ ~ \aPlJ{heYi)f ~ 14.3 MW/cm^. If the inci- 
dent pulse intensity is well below this limiting one the 
density of excitons in GaAs waveguide is small enough 
and both Eq. ^ and the developed theory as a whole 
remain valid. It should be noted here that the dam- 
age threshold for the chosen material is even higher than 
-^iim- about 40 MW/cm^ for a 20 ns incident pulse.l^ 
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As this threshold grows rapidly with the reduction of in- 
cident pulse duratiorP^ one can expect that it is much 
higher for ps-pulses considered in this work. Knowing 
the nonlinear coefficient a one can also easily reproduce 
the input pulse intensities, which are necessary to form 
solitons in the excitonic waveguide. According to the 
used scaling of electric field (Eq = 1 / {^y atQatQEB)) , for 
example, for solitons with a temporal width to = 2/a (see 
Fig.p| these intensities are: ^10— 100 kW/cm^ for inci- 
dent pulses spectrally centered inside the polaritonic gap 
and ~ 0.1 — 1 kW/cm^ for those spectrally centered out- 
side the gap (it is clear that in both cases these intensities 
are well below the limiting ones). 

Nonlinear effects such as solitons formation could be 
observed, for example, in silica optical fibers possessing 
a Kerr nonlinearity only at much higher pulse intensities^ 
(< 1 MW/cm^). The source of this large difference in in- 
tensities can be found by the proper comparison of the 
nonlinear coefficient a at the cubic term in Eq. ([2| with 
the corresponding coefficient of the Kerr effect in silica 
fiber^ (which is actually the nonlinear refractive index 



,Si02 
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m /W). As we will see below the exci- 



tonic nonlinearity is much stronger than the Kerr one, 
however it has a resonant nature - it decreases rapidly 
with the increase of detuning. This can be demonstrated 
by the example of plane wave propagation. In this case, 
according to Eq. ^ the amplitudes of polarization and 
electric field are related by the following expansion up to 
the third order: 



MF 



E - acsE , Q!off 



(Q£BrMF)^ri 



This dependence clearly shows that the factor multiply- 
ing the cubic nonlinear term decreases rapidly when the 
detuning Awx grows. In other words, if the incident pulse 
is spectrally centered far enough from the exciton res- 
onance the considered nonlinearity becomes negligible, 
and the background Kerr nonlinearity of the medium 
starts to play the main role. However, near the reso- 
nance (Awx ~ 10a) for reasonabl^^H low-temperature val- 
ues of /i7x — 0.03 — 0.3 meV one can estimate Re a^e ^ 
10^^ — lO^^nj'*^^, that clearly indicates the domination 
of excitonic nonlinearity in the vicinity of resonance. 

It is also important to note that excitonic media effi- 
ciently convert light into excitons, thus being optically 
highly absorptive near the exciton resonance. Therefore, 
the incident pulse spectrum central frequency should not 
be too close to Wx. One can estimate the corresponding 
absorption length a^^H 



I = 



c Re ^/e 
uj Im e ' 



Aw^ 



MFi FM 

- j7x 



(11) 



where e is the polariton dielectric function in the linear 
regime. In the case of GaAs, for reasonable values of the 
detuning (a few times a), this length is of the order of 
several microns. This agrees well with the micron-scale 
lengths, which are completely sufficient for the revela- 
tion of nonlinear optical effects in excitonic waveguides 




FIG. 5. Schematic drawing of a sandwiched semiconducting 
waveguide. The realistic length parameters are indicated for 
the case of GaAs (central layer) surrounded by AlGaAs. 



(see Sec. III). However, waveguide lengths of just a few 



microns are not sufficient to actually guide light, and it 
is unlikely that one can produce such a short waveguide 
made of pure semiconducting material. One possible so- 
lution is to use a composite waveguide made of material 
without excitonic resonances in its optical response (e.g., 
silica) doped with macroscopic domains of semiconduc- 
tor (for example, nanoparticles) . But then a different 
theory based on the Maxwell-Garnett formalism should 
be developed. Another way is to use a waveguide in the 
form of sandwiched semiconducting heterostructure with 
central working layer made of semiconductor with an ex- 
citonic feature surrounded by the semiconducting layers 
made of a different material (see Fig. [S]). This cladding 
material should have no excitonic resonances in the same 
range of frequencies than that in the working layer so 
that the propagating light will not be affected by the ex- 
citonic nonlinearity, but should possess almost the same 
refractive index so that the optical guided mode is well 
spread outside the working layer, and the photon-exciton 
coupling is reduced accordingly. 

Because of this special structure of the proposed 
waveguide, the modal distribution F for the electric field 
will be much wider than the distribution M for the po- 
larization field. Therefore, the transverse integral Ffm in 
the main system ([3| becomes much smaller than unity, 
while still Fmf ^ 1- For example, by using the realis- 
tic AlGaAs/GaAs sandwiched waveguide transverse pa- 
rameters, which are indicated at Fig. [5j one can obtain 
Ffm ^ 1 and then get from Eq. ( [ll| ) for Awx ^ 10a an 
absorption length I ^ 100 /im. Of course, as it follows 
from system ([3| such a small Ffm leads to a weak cou- 
pling between the electric and polarization fields and thus 
to a weak influence of the excitonic nonlinearity on the 
electric field, but this is compensated by the long propa- 
gation distance. Also the reduction of Ffm results in de- 



crease of the gap width Ap (see Sec. IV), therefore longer 



incident pulses should be used to have a sufficiently nar- 
row spectrum and efficiently deliver the pulse energy into 
the polaritonic gap. 

As a last point we should also note that low- 
dimensional semiconductor structures, such as semicon- 
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ductor microcavitiesp' certainly have a number of ad- 
vantages with respect to bulk semiconductors considered 
here (e.g., larger exciton binding frequency, which trans- 
lates into larger optical nonlinearities). However, waveg- 
uides based on a bulk semiconductor are preferable for 
observation of nonlinear processes associated with soli- 
tons because they are weakly affected by diffraction, 
which is hard to overcome in microcavities,'^^ and also 
because such waveguides can be more easily produced 
and tested than guiding structures in microcavities. A 
possibility to properly combine the excellent waveguid- 
ing properties of an optical fiber with optical features 
of pure bulk semiconductors in a single structure has 
been very recently demonstrated - a silica optical fiber 
with a semiconducting ZnSe core has been fabricated.!^ 
This should also stimulate further experimental investi- 
gations of nonlinear optical properties of bulk semicon- 
ductor waveguides. 

VI. CONCLUSIONS 

In conclusion, we have carried out a MI analysis of con- 
tinuous waves propagation in semiconducting excitonic 
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